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Abstract 

Decay or annihilation products of superheavy dark matter (SHDM) could be 
responsible for the end of the Ultra-High Energy Cosmic Ray (UHECR) spectrum. 
In this case, the south array of the Pierre Auger Observatory should observe in 
the future a significant anisotropy of UHECR arrival directions towards the galac- 
tic center. Here we use the already existing data of the SUGAR array to test 
this possibility. If decaying SHDM is distributed according a Navarro-Frenk- White 
(NFW) dark matter profile with core radius R c = 15 kpc and is responsible only for 
UHECRs above ~ 6 x 10 19 eV, i.e. the AGASA excess, then the arrival directions 
measured by the SUGAR array have a probability of ~ 10% to be consistent with 
this model. By contrast, the model of annihilating SHDM is disfavoured at least at 
99% CL by the SUGAR data, if the smooth component of the DM dominates the 
signal. 

PACS numbers: 98.70. Sa, 14.80.-j 



1 Introduction 

Protons accelerated by uniformly distributed extragalactic astrophysical sources would 
be a perfect minimal explanation of the UHECR data above 10 19 eV. However, protons 
with energy E > 4 x 10 19 eV loose quickly energy due to pion production on cosmic 
microwave background photons. Thus the proton spectrum should show the so-called 
Greisen-Zatsepin-Kuzmin (GZK) cutoff pQ, which is not observed by the AGASA ex- 
periment f2J. Moreover, if the small-scale clusters in the arrival directions of UHECR 
measured by AGASA j2| are due to point-like sources, one can estimate their num- 
ber 0. This number is so small that the nearest source should be located at the dis- 
tance -R m in ~ 100 Mpc ^6J. This means that the GZK cutoff is exponentially sharp at 
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E rj 6 x 10 19 eV and even the data of the HiRes experiment [7] are inconsistent with 
the expected proton spectrum jHj. This inconsistency becomes even stronger if BL Lacs, 
which show a statistically significant correlations with the arrival directions of UHECR 
with energy £~4-6x 10 19 eV jS], are sources of UHECR. 

A possible solution to this problem would be the existence of superheavy dark matter 
(SHDM) jHllini- Superheavy particles with mass M x ~ 10 13-14 GeV can be naturally 
produced during inflation and would be today the dominant component of dark mat- 
ter [TT]. Such particles will concentrate in galactic halos and the secondaries from their 
decay could be responsible for the highest energy cosmic rays. It has also been suggested 
that not decays but annihilations of SHDM particles produce the observed UHECRs [T7] . 

By construction, this model has two clean signatures: the dominance of photons at 
the highest energies |9| and an anisotropy of the arrival directions with an increased flux 
from the Galactic center [T2*l 113 j . Unfortunately, both signatures are not very clean for 
the present experiments. Indeed, at 95% C.L., ~ 30% of the UHECR above E > 10 19 eV 
can be photons |14j . which means that still most of UHECRs with E > 4 x 10 19 eV can 
be photons without any contradiction to the experimental data. Since experiments in 
the northern hemisphere do not see the Galactic center, they are not very sensitive to a 
possible anisotropy of arrival directions of UHECR from SHDM. In contrast, the Galactic 
center was visible for the old Australian SUGAR experiment. 

The anisotropy of the arrival directions using data from the full sky was discussed in 
Refs. [T5| ITE] . Reference compared the flux from the Galactic center to the one from 
the anti-center and found them to be comparable. Similarly, the full-sky harmonic anal- 
ysis including AGASA and SUGAR data from Ref. ^H] found no significant anisotropy. 
In this work, we use a two-component energy spectrum of UHECRs consisting of pro- 
tons from uniformly distributed, astrophysical sources and the fragmentation products of 
SHDM calculated in SUSY-QCD. We compare their expected arrival direction distribu- 
tion to the data of the SUGAR experiment using a Kolmogorov-Smirnov test. Contrary 
to the harmonic analysis, this test allows to quantify directly the (dis-) agreement of the 
measured distribution of arrival direction with the expected one in the SHDM model. We 
consider both decays and annihilations of SHDM. 

The paper is organized as follows: In Section II we discuss the status of the SUGAR 
data. The contribution of SHDM to the UHECR spectrum is discussed in Section III. 
In Section IV we perform an harmonic analysis of the arrival directions measured by the 
SUGAR experiment. Then we use a Kolmogorov-Smirnov test in one and two dimensions 
to check the consistency of the SUGAR data with the probability distribution of arrival 
directions expected for SHDM in Section V. Finally we conclude in Section VI. 

2 Assessment of the SUGAR data 

In order to use the SUGAR energy spectrum correctly we compare their data given in 
Ref. [20J to the energy spectrum measured by the AGASA j2] and HiRes experiments [7] 
in Fig. [TJ Rescaling the SUGAR energies calculated with the Hillas prescription by 15% 
downwards, -EWias —* ^HiUas/1-15, makes their data consistent with the ones from AGASA. 
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Figure 1: UHECR spectrum measured by the AGASA, HiRes and SUGAR experiments. 
We scale the SUGAR spectrum down by E/1.15 and the HiRes spectrum up by 1.25 x E 
using the AGASA spectrum as reference. The overall normalization of the spectra has 
only a weak influence on our results. 



The same is true for the HiRes spectrum if the energy is rescaled up-wards by 25% 1 
We have chosen arbitrary the AGASA spectrum as reference, but changing the overall 
normalization of the spectra affects our results only weakly. As seen from Fig. ^ the 
SUGAR spectrum has the ankle at the correct place around E « 10 19 eV and is consistent 
with the AGASA spectrum in the whole energy range. In particular, the SUGAR spectrum 
also does not show the GZK cutoff at the highest energies. 

The rescaling of the SUGAR data downwards by 15% should be compared to the recent 
reevaluation of the energy conversion formula used in the Haverah Park experiment [22] • 
In this reference, the relation between p(600) and the primary energy has been recalculated 
using QGSJET [21] and compared to the original relation suggested by Hillas. The new 
calibration results in ~ 30% lower primary energies. 

The SUGAR experiment was a very sparse array of detectors and its energy deter- 
mination of each single event was therefore rather unprecise. Thus, we shall use as a 
statistical test later on a method which relies only on the total flux measured by the 
SUGAR array, but uses not the energy of each single event. Since after the rescaling of 
the energies measured by SUGAR, -Emilias — > -^Hiiias/1-15, its measured flux is consistent 
with newer experiments like AGASA and HiRes, we conclude that on average the energy 
determination in the SUGAR experiment was reliable. 

1 The aperture of HiRes is energy dependent; the rescaling we perform should be seen therefore just 
as a crude approximation. 
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The energy conversion formula used in SUGAR to connect the measured muon number 
with the primary energy assumes that the primary is a hadron. For photon primaries, 
predicted to be dominant in the SHDM model, the muon content of the shower is smaller 
by a factor 5-10 [23 . Thus the energies of photon events is expected to be underestimated 
by the SUGAR experiment. The SUGAR spectrum shown in Fig. 1 would be unchanged 
at energies E < 5 x 10 19 eV, i.e at energies where all three experiments agree after 
rescaling. 

The angular acceptance r](5) as function of declination 5 averaged over time of an 
experiment at geographical latitude b {b = —30.5° for SUGAR) observing showers with 
maximal zenith angle 6 max is 



r](S) cx / da cos(#) oc [cos(6) cos(5) sin(a max ) + a max sin(6) sin(<5)] (1) 
Jo 

where 

t _ cos(6» max ) - sin(6) sin(5) 



cos(6) cos(<5) 



and 




for - 1 < f < 1 , 

"max = { 71 for f < -1 , (3) 

for £ > 1 . 

We have checked that the zenith angle distribution of the SUGAR events agrees with 
the theoretical predicted one, d/Vth oc d6sin(6) cos(#), above E > 4 x 10 19 eV. At lower 
energies, the acceptance of the experiment becomes energy dependent and deviations from 
dN t h start to grow. 



3 Superheavy dark matter contribution to UHECR 
spectrum 

We fix the contribution of SHDM to the total UHECR flux following the assumptions of 
Ref. [18J: we assume that no galactic astrophysical sources contribute to the cosmic ray 
flux above 10 19 eV and that the extragalactic cosmic ray flux can be characterized by 
an injection spectra of protons with a single power law, j ex {E) oc E~ a . For the choice 
of a = 2.7, this energy spectra modified by redshift, e + e~ and pion production fits very 
well the measured spectra below E < 6 — 8 x 10 18 eV ^B]- The only difference with ^H] 
is that we take into account that total number of sources is small |6,, if the small scale 
clusters measured by the AGASA experiment are due to point-like sources. The AGASA 
data favor as minimal distance to the nearest source -D min ~ 100 Mpc 6J. Therefore, the 
contribution of protons from extragalactic sources has a sharp cutoff. For the calculations 
of the proton spectrum we used the code [2E]- We use then the SUSY QCD fragmentation 
functions D(x, M x ) of superheavy particles with mass Mx calculated in Ref. JH] to model 
the flux j'dm(-S) oc D 7 (x, M x )- The total UHECR flux is thus 

j{E) = (l-e)j ex {E) + ejv U {E). (4) 



4 




10 iu Iff" 1U 1U 

E [eV] E m 



Figure 2: UHECR spectrum measured by AGASA experiment. Protons from extragalactic 
sources contribute below the GZK cutoff, photons from SHDM decays mainly above the 
GZK cutoff; for hard, 1/E 2 - 3 , (left panel) and soft, 1/E 2 - 7 (right panel) injection spectrum 
of extragalactic protons. 



We fix the constant e determining the relative contribution of SHDM to the UHECR 
flux by a fit of j(E) to the AGASA data j2J. In Fig. |2J we show our fits for the case 
of a harder 1/E 2 ' 3 (left) and a softer 1/E 2 ' 7 (right) injection spectrum of extragalactic 
protons. In the first case, the contribution of SHDM to the UHECR spectrum below the 
GZK cutoff is minimal and starts to dominate only at highest energies E > 6 x 10 19 eV. 
For this choice of injection spectrum, the contribution of galactic sources dominate for 
E < 10 19 eV. In the second case, SHDM gives a larger contribution at lower energies 
E < 6 x 10 19 than before, and again starts to dominate at E > 6 x 10 19 eV. Note, that 
most UHECRs above E > 6 x 10 19 eV should be photons in this model, but this does not 
contradict the rather weak existing bound of 30% of photons at E > 10 19 eV |14j . 

In the following, we shall use conservatively the case of the harder 1/E 2 - 3 injection 
spectrum if not otherwise stated. Then the contribution of SHDM to the UHECR spec- 
trum is fixed by Eq. (£Q). 



4 Harmonic analysis 

In order to compare SUGAR data with an uniform distribution typical for extragalactic 
sources we have performed an one-dimensional harmonic analysis. As usual we sum 

2 n 2 n 

dk = — cos (k(fi) and bk = — sin(/c0) (5) 

n o=l n 0=1 

over the n data points. 

The amplitude and and phase <fik of the k.th harmonic are given by 

r k = \l a \ + b 2 k and 4> k = arctan(6 fc /a A .) (6) 
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Table 1: Direction to the signal in right ascension and chance probability of the k.th 
harmonics to arise from an isotropic distribution; for different cuts in energy E and zenith 
angle 9 < 55°. 

with chance probability 

p ch = exp (-nrj|/4) . (7) 

The direction to the signal is <fi = k(f>k- 

Results of a harmonic analysis in right ascension a depending on E m[n and # max are 
given in Table ^ The results for all harmonics show generally good agreement with 
an isotropic distribution for any cutoff energy we have used. Only the third harmonics 
shows some anisotropy, in particular at the highest energies, E > 8 x 10 19 eV; however 
its phase does not points towards the galactic center (lying at a = 266°). Generally, 
all harmonics point instead towards a ~ 130°. Reference derived the probability 
distribution (pdf) that a data set with phase <pi and amplitude r± is drawn from an 
arbitrary pdf. However, we are not aware of a generalization to higher harmonics and, in 
particular, of a method to combine the information content of several harmonics. In the 
next section, we use therefore a Kolmogorov-Smirnov test to quantify the (dis-) agreement 
between the expected distribution of arrival direction and the SUGAR data. 



5 Kolmogorov-Smirnov tests 

The pdf to detect an event with energy E and arrival direction a, 5 is a combination of 
the isotropic extragalactic and the SHDM flux, 



p(E, a, 5) oc r)(d) 



Jcx 

(E) +j DM (E) 



ds nv M (r(a,6)) 



where s max = Re cos 6 + \J R\ — R% sin 2 6 is given by the extension R h = 100 kpc of the 
DM halo and 9 is the angle relative to the direction to the GC. As explained, the relative 
size of the two contributions is fixed by the fit to Eq. (HJ). 

For the two-dimensional test, we have integrated Eq. (JHJ) over energy, 



P2d(ac, S) 



dE p(E, a, 5) 



(9) 
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where E m - m and E max are the minimal and maximal energy considered in the UHECR 
spectrum. We have used £ max = io 21-22 eV, but the results do depend on very weakly 
on the exact value. By contrast, the value of E m[n has a strong influence on the results 
obtained. 

For the one- dimensional test, we have integrated Eq. Q over the declination, 

Pcn{a) = r 2 d5 cos 5 P 2d {a,5) . (10) 

J-w/2 

In the standard one-dimensional Kolmogorov-Smirnov (KS) test, the maximal differ- 
ence D between the cumulative probability distribution function P(x) = J dx'p(x') and 
the cumulative distribution function of the data, 

S(x) = -^29(xi-x), (11) 

is used as estimator for the belief that the data are drawn from the distribution p. A 
variant of this test which is equally sensitive on differences for all x and is especially well 
suited for data on S 1 uses instead of D the symmetric estimator 

V = D + + D_ = max[S(x) - P(x)\ + max[P(x) - S(x)\ . (12) 

The significance of a certain value of V is calculated with the formula given in Since 
the exposure of a ground-array experiment is uniform in right ascension a, we use a as 
variable in the one-dimensional KS test. More exactly, we use as pdf Eq. (JHJ) integrated 
over dE and d5 cos 6. 



0.1 



0.01 



0.001 



1 o 




55 degrees. _ — 




0.8 


2c 




0.6 


/ 

/ 

45 degrees 




0.4 




3a 


0.2 








10 



15 20 25 

R c [kpc] 



30 




Figure 3: Left: Consistency level of the SUGAR data with SHDM distributed according 
a NFW profile as function of the core radius R c ; SHDM is assumed to be the source 
of all UHECRs above E min = 8 x 10 19 eV. Right: Comparison of S(a) and P{a) for 
R s = 15 kpc, E n 



. mul - 8 x 10 19 eV and # max = 45°. 



As simplest test, we assume that all SUGAR events above E min are produced by 
SHDM. Thus we compare S(a) with P(a) = Jq da'PcR^a'). The result is shown for 



7 



E > 8 x 10 19 eV in Fig.|SK for two different values of the maximal zenith angle, # max = 45° 
and 6* max = 55°. While for 6* max = 45° SHDM is disfavoured at the two sigma level for 
realistic values of the core radius, R c ~ 20 kpc, the SUGAR data have for the choice 
of #max — 55° a rather large probability p to be consistent with the SHDM hypothesis, 
p ~ 20%. In Fig. Eb; we compare the two commulative distributions S(a) and P{a) 
for R s = 15 kpc, E min = 8 x 10 19 eV and 9 max = 45°. Inspecting S(a) makes it clear 
that the data in this case are not uniformly distributed but clustered around a ~ 130° 
and a ~ 350°. Since none of these two directions coincide with the position of the GC, 
this data set disfavours the SHDM hypothesis more strongly than one would expect for 
uniformly distributed events from extragalactic sources. However, one should use a rather 
low value of E m \ n to minimize the uncertainties in the SUGAR energy determination and 
we will therefore not rely on these results. 

We consider therefore next as more realistic test the case that both SHDM and ex- 
tragalactic sources contribute to the UHECR spectrum. Then the dependence of p on 
E m \ n should be diminished. More exactly, one would expect in the case that the SHDM 
hypothesis is disfavoured by the data that decreasing E m i n first decreases p. This de- 
crease of p should continue down until E ~ (3 — 4) x 10 19 eV, i.e. until a point where the 
signal-to-background ratio becomes considerably smaller than one. Decreasing E min even 
further should result in a increase of p because now practically all new events are from 
extragalactic sources. 
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Figure 4: Left: Dependence of the probability on the energy cutoff in the SUGAR data for 
decaying SHDM. Right: Two-dimensional KS test give results similar to one-dimensional 
test. 



In Fig. Hk, we show the dependence of the probability on the energy cutoff for R c = 
15 kpc and # max = 55°. The two thick solid lines show p for a combination of SHDM and 
uniform sources according Eq. the upper one corresponds to an injection spectrum 
1/E 2 ' 3 , the lower one to an injection spectrum 1/E 2 ' 7 . The behaviour of p suggest that 
the minimum for E m i n ~ 7 x 10 19 eV is a fluctuation similar to the maximum around 
-Emin ~ 5 x 10 19 eV. In the range E min ~ (3 — 4) x 10 19 eV, the fluctuations adding an 
additional event are relatively small. Therefore, we consider the probability in this range 
as more reliable indicator for the consistency of SHDM with the SUGAR arrival directions; 
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we conclude that the SUGAR data have the probability p = 5 — 20% to be consistent 
with the SHDM depending on the injection spectrum of the extragalactic protons. 

The thin solid line shows how consistent the SUGAR data are with an isotropic dis- 
tribution. This distribution has also a minimum around £ min ~ 7 x 10 19 eV where 
the events cluster around two arrival directions. After including more low-energy data, 
the SUGAR arrival directions are consistent with an isotropic distribution. Finally, the 
dashed line shows the consistency of the SUGAR data with the assumption that all 
UHECR events above E min are from SHDM. It is clear that only values of E min above 
E m i n ~ 4 x 10 19 eV are compatible with the SUGAR data. Similar, the spectral shape of 
the flux in the SHDM model allows a dominance of SHDM in the UHECR spectrum only 
above E > 6 x 10 19 eV |19j . In Fig. ^3 we compare results from one- and two-dimensional 
KS tests (of a and 5) as function of the energy cutoff and find that they give rather similar 
results. 
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Figure 5: Dependence of the probability on the energy cutoff in the SUGAR data for 
annihilations of SHDM. Left: for core radius R c = 15 kpc and different e determining the 
SHDM contribution. Right: for several core radii R c ; assumes that all events above E are 
from SHDM. 



Finally, we consider the model where not decays but annihilations of SHDM particles 
produce the observed UHECRs 17 j. In the original version of this model it was suggested 
that the flux of the clumpy component dominates over the one from the smooth SHDM 
profile by 3 orders of magnitude. On the other side, it was shown in a recent paper [3U| 
that the contribution of the clumpy component can be just a factor few larger than the 
one of the smooth component. Moreover, the newest numerical calculations show that 
the contribution of clumps is even subdominant and that it is very unlikely that a nearby 
clump will outshine the Galactic center |3T] . Because of the arguments above, we assume 
that the clumpy part of the SHDM gives a subdominant contribution to the UHECR 
flux. In the opposite case our results for the SHDM model with annihilations will be less 
significant, depending on the relative contribution of the two components. 

Since the flux is now oc Hq M , the anisotropy in this model is much stronger than for 
decaying SHDM. This can be clearly seen in Fig. where we show the dependence of 
the probability of annihilating SHDM on the core radius R c assuming that all events 
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above E are from SHDM. Even for core radii as large as 30 kpc, annihilating SHDM 
is disfavoured by two sigma. Figure shows similar to Fig. H^l the dependence of the 
probability on the energy cutoff for R c = 15 kpc and # max = 55°. The two thick solid lines 
show p for a combination of SHDM and uniform sources according Eq. (@J); the upper 
one corresponds to an injection spectrum 1/E 2 ' 3 , the lower one to an injection spectrum 
1/E 2 ' 7 . Depending on the injection spectrum of extragalactic protons, annihilations of 
SHDM are disfavoured by the SUGAR data between 3 and 4a. 



6 Conclusions 

In this paper we have tested the consistency of the SHDM model with the SUGAR data. 
In order to use the SUGAR data, we have compared its energy spectrum to the one of 
AGASA and found that they are compatible after rescaling down the SUGAR energies by 
15%. We have assumed that the energy spectrum in the region 10 19 eV < E < 6 x 10 19 , 
i.e. between ankle and GZK cutoff, is dominated by protons coming from uniformly 
distributed extragalactic sources. After fitting the relative contributions of SHDM decay 
products and extragalactic protons to the AGASA data, we have performed Kolmogorov- 
Smirnov tests of the SUGAR data. As result we have found that SUGAR data are able 
to disfavour strongly extreme case like annihilations of SHDM without clumps (5a) or 
decaying SHDM (99% CL) assuming their contribution to the UHECR flux dominates 
down to E = 4 x 10 19 eV. The phenomenologically most interesting case, decaying SHDM 
dominating the UHECR spectrum only above E > 6 x 10 19 eV, is consistent with the 
SUGAR date with 5-20% probability. Thus the SUGAR data do not disfavour strongly 
this model but they neither support it. A statistically significant test of this model can 
be done by the Pierre Auger Observatory. 
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